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Abstract 

A process of single diffractive dissociation off nuclei is considered on a basis of solutions to 
the nonlinear evolution equation. The relevant saturation scales Q^Ax, xq) are determined 
and their dependences on Bjorken x, atomic number A, and minimal rapidity gap Yq = 
In 1/xq are investigated. The solutions are shown to possess a geometrical scaling in a broad 
kinematic region for x well below xq. 

The ratio a diff /a tot is computed for several nuclei. We predict that at x ~ 10 4 this 
ratio is of the order 25% which is much larger compared to the one of proton. This result 
indicates a possibility to observe a very strong nuclear shadowing. 
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1 Introduction 



During the last years there has been a significant growth in the interest to a new phase of QCD 
associated with high parton density 0,0. This interest is mostly related to availability of new low 
x e — p DIS data. Another source of information on QCD dynamics at high parton density is due 
to nuclei which can provide high density effects at comparatively lower energies. The researches 
on nuclear shadowing has been recently accelerated due to the start of the RHIC collider. 

In the present paper we concentrate on a process of single diffraction dissociation off nuclei which 
plays an important role in revealing QCD dynamics at high parton density. Diffractive inclusive 
production in DIS is believed to be very sensitive to shadowing effects PJ and, in fact, is a measure 
of these effects (due to the AGK cutting rules ||). We are going to investigate the cross section 
of diffractive dissociation at fixed impact parameter and as a final result compute the ratio of the 
total inclusive diffraction to the total inclusive cross section in DIS off nuclei. 

The total deep inelastic cross section is related to the dipole cross section 

a{x,Q 2 ) = J d 2 r ± Jdz\^*{Q 2 ;r ± ,z)\ 2 a dipolc {r ± ,x), (1.1) 

where the QED wave functions \1/ 7 * of the virtual photon are well known [EL BL PJ. The dipole 
cross section is given by the integral over the impact parameter b: 



Cdi P oie(?~±,:r) 



d 2 b N(r ± ,x;b) 



;i.2) 



where N stands for imaginary part of the dipole-target interaction amplitude for dipole of the 
size r_i_ scattered elastically at the impact parameter b. This function is a subject to a nonlinear 
quantum evolution equation first derived by Balitsky and Kovchegov (BK) ||, |[| . The BK equation 
was studied both asymptotically @ and numerically |L0|, |TTJ. For the purposes of the present 



paper we will use a numerical solution of this equation obtained in Ref. \V1 
The total cross section of single diffractive dissociation is similarly defined: 



a diff (x,x ,Q 2 ) = J d 2 r ± J dz (Q 2 ;r ± , z)\ 2 a^(r ± ,x,x ), (1.3) 

and 



a 



diff 
dipole 



(r±,x,x ) 



drb N u (r x ,x,x ;b) 



;i-4) 



The function N D is a partial cross section for the dipole-target diffractive scattering with the 
minimal rapidity gap Y = Ib.1/xq. A nonlinear quantum evolution equation for N D was derived 
in Ref. [R| and then rederived in Ref. [14]: 
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x /„ rf 2 x 2 ^-[ 2 iV D (x 02 , y, >o; b - |x 12 ) + iV D (x 02 , y, Y ; b - §x 12 )iV D (x 12 , y, Y ; h - |x 02 ) 

-4 iV D (x 02 , y, F ; b - §x 12 )iV(x 12 , y; b - |x 02 ) + 2 iV(x 02 , y; b - |x 12 )iV(x 12) y;b- |x 02 )] . 

(1.5) 

The lhs. of Eq. (|1.5|) is a partial cross section of the diffractive dissociation for dipole of the size 
r± = x i = |x — xi| and rapidity Y = lnl/x. The rhs. of Eq. ( |1.5| ) describes quantum evolution 
in which the original dipole first splits to two dipoles and then the latter scatter off the target. 
The ultraviolet cutoff p is defined to regularize the integral, but it does not appear in physical 
observables. The following numerical results are checked to be indpendent on a choice of p. 

When the rapidity fills the whole rapidity gap (Y = Y ), diffraction is reduced to pure elastic 
interaction 



N D (r ± ,x ,x ;b) = N 2 {r ± ,x ;b) . 



;i.6) 



Eq. (|1.6| ) serves as initial condition to the evolution ( |i~5] 



For a proton target a numerical solution of Eq. ( |1.5|) was found and investigated for the first 
In the following paper the ratio a diff /a tot was computed and in a certain 

In 



time in Ref. [16 
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kinematic domain shown to be energy independent in agreement with the HERA data 
the present paper we report on the numerical solution of Eq. ( |1.5| ) for nucleus targets and repeat 
the program of Refs. [16, 17]. 

The paper is organized as follows. In the next Section (2) we present solutions of Eq. ( |1.5| ) for 
several nuclei. Section 3 is devoted to the determination of the saturation scale and its properties. 
Geometrical scaling is studied in Section 4. The ratio a dl ^ /a tot is computed in Section 5. Section 
6 present a theoretical discussion of the saturation scale. The last Section (7) is concluding. 



2 Solution of the non-linear equation 



In this section we report on the numerical solution of Eq. ( |1 . 5| ) for several nuclear targets. We 
consider six real nuclei: Ne2o, C«4o, Zn 70 , Moioo, A^iso, and Auigj. All the details about nuclear 
profile functions are borrowed from Ref. [T^| and are summarized in our paper ||12|| . As in a series 



of our previous papers, solutions to Eq. ( |1.5|) are obtained by the method of iterations proposed 
in Ref. [[□]]. The constant value for the strong coupling constant as = 0.25 is always used. For 
the function N, which is an input in Eq. (|1.5| ), we use a solution of the BK equation obtained 



in Ref. [12||. The solutions are computed for 10 < xq < 10 and within the kinematic region 



10 ' < x < Xq and transverse distances up to a few fermi. 
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The function iV^Q is formally a function of four variables: the energy gap x$, the energy variable 
x, the transverse distance r±, and the impact parameter b. In order to simplify the problem we 
will proceed similarly to the treatment of the 6-dependence of the function Na [Q an d N D for 
proton |T3 . Namely, we use the ansatz which preserves the very same 6-dependence as introduced 
in the initial conditions: 



N°(r ± ,x,x ;b) = (1 - e -*2(^o.rj^W/5 A (o) )2j 



(2.7) 



with 



K A {x,x ,r± 



) = -ln(l - >jN°{r ±1 x,x )). 



(2.1 



The function (r±, x, Xq) is a solution of Eq. ( |1.5| ) but with no dependence on the forth variable. 
The initial conditions for N^(r±, x, Xq) are set at b = 0. 

Fig. [I] displays the solution as a function of the transverse distance for several values of x 
and at fixed xq = 10~ 2 . The obtained curves reproduce the very same pattern as in the case of a 
proton target . The results for various nuclei can be used in order to study the A dependence 
of the diffraction dissociation. In agreement with all expectations, the unitarity bound = 1 
is reached first by the most heavy nucleus. The dependence of the solutions on the gap variable 
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Figure 1: The function N®{xq = 10 2 ) is plotted versus transverse distance. The curves correspond to 
various nuclei (from up to down): Au, Nd, Mo, Zn, Ca, Ne. 



Xq is very weak and quite similar to the proton case of Ref. [16 . 



3 Saturation Scale 



Determination of the diffractive saturation scale (A, x, xq) from the solution N% is a subject 



of this section. Following the same spirit of our previous works [11 , V2\ 20 , [Tj| we introduce 



several definitions of the saturation scale while a variety of thus obtained results will indicate an 
uncertainty of the determined values. 

For the step like function it is natural to define the saturation scale as a position where = 1/2: 
lr The subscript A indicates a relation to a nuclear target with atomic number A. 



3 



Definition (a): 



N%(R?,x,x ) = 1/2, Q u s = 2/R 



D 



(3.9) 



The equality between the saturation radius R® and the saturation momentum is motivated 
by the double logarithmic approximation. Though this approximation is formally not justified, 
we still believe it to make reliable estimates provided Qf is large enough. The definition ( ^.9| ) 
is analogous to the one proposed in Ref. 0, N{2/Q s ,x) = 1/2. If we recall that N% = N\ at 
x = xq and postulate Qf (A, x , xq) = Q s (A,xq) then we should require 



Definition (b): 



N%(2/Q»,x,x ) = 1/4. 



(3.10) 



An alternative definition of the saturation scale could be one motivated by the Glauber-Mueller 
formula || for Na(t±,x) 



N A (r±,x,x ;b) 



with 



K A 



3 7r 2 as 



V A xGiyX , — 2~ 



r 



4 

2~ I ' 



(3.11) 
(3.12) 



and xG standing for the gluon density of a nucleon. 
• Definition (c): 



k%(2/Q?,x,x ) = 1/2 



(3.13) 



The obtained saturation scales are depicted in Fig. |2] for xo = 10 2 and in Fig. for xq = 10 3 . 
The saturation scales are practically Xq independent. The definition (d) is related to scaling 
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Figure 2: The saturation scale Q®{xq = 10 2 ) is plotted versus x. The curves correspond to (from up 
to down) Au, Mo, Ca, and Ne. 

properties of the function and will be discussed in the next section. 

It is important to learn about x and A dependencies of the saturation scale. To this goal, we 
assume the following parametrization: 

Qf(A, x, x ) = Qf (x ) A p x- x . (3.14) 
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Figure 3: The saturation scale Q®(xq = 10 3 ) is plotted versus x. The curves correspond to (from up 
to down) Au, Mo, Ca, and Ne. 



Nuclei \ x 


io- 7 


10~ 6 


10~ 5 


IO" 4 
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Light 


0.15 


0.20 


0.24 


0.29 


0.35 


Heavy 


0.15 


0.19 


0.22 


0.25 


0.29 


All 


0.15 


0.19 


0.23 


0.27 


0.32 



Table 1: The power p(x) for various val- 
ues of x. 



In fact, the parametrization ( 3.14 ) is a good approximation for the obtained values of the satu- 
ration scales with 

A = 0.37 ± 0.04, 

and p given in Table [TJ. All the values in the table are central values given with errors less than 
10%. 

The Xo-dependence of the saturation scale is very weak. If we try to use the power law parametriza- 



tion Q 



D 



ft 



then (3 = 0.055 ± 0.035 . The large error reflects a significant x-dependence of (3 



as well as its sensitivity to a saturation scale definition. 

It is important to stress that the obtained power A coincides with the corresponding power of the 
saturation scales Q s [|TJ] and Qf [16). The error obtained for A is purely numerical. We have 
to admit that two independent collaborations |1C| found A ~ 2 as which exceeds our value by 
about 25%. This discrepancy may result from different definitions of the saturation scale. We 
remind that in our approach we actually compute the saturation radius and not a momentum. A 
transition to the latter is likely to be more complicated compared to the relation Q s = 2/R s . Yet, 
our estimates of the saturation scale based on the solutions in momentum space do not indicate 
any dramatic change in A. In addition, we have want to emphisize that we use completely diffrent 
initial conditions compared to Ref. 



It was shown in Ref. 



that in fact solutions to the 



non-linear evolution equation crucially depend on a choice of initial conditions. 



4 Geometrical Scaling 



In Ref. [20] the function N was shown to display a phenomenon of geometrical scaling while Ref. 



16 presents a detailed analysis of the scaling for diffractive dissociation off proton. As can be 
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expected, the function displays the very same property and in this section we give a brief 
illustration of the phenomena. In the saturation region the scaling implies the amplitude to be a 
function of only one variable r = r\ - (Qf(A, x, x )) 2 : 

N°(r ± ,x,x ) = N°(t) (4.15) 

Let us define the following derivative functions while the second equalities hold if the scaling 
behavior ( 4.15|) is assumed: 



, T Dr „ \ 9N% dN% <91n(Qf) 2 ft „. 
N°(r ± ,A,x,x ) = --g± = -+r-£d-, (4.16) 

N?(r x ,A,x,x o )=ri?0- = ^T, (4.17) 

a \ dN$ dN% <9hi(Q?) 2 fA ,. 
$l(r ± ,A,x,x ) = -— f = — A r Ws . 4.18 
oYq dr a In x 

If the scaling behavior (|4.15|) takes place indeed, then both ratios /N® and $t/N® are r_|_ 
independent functions. 

Let us first consider the scaling with respect to x. Fig. |4] presents the derivatives and as 
functions of transverse distance r± at fixed Xq = 10 -2 . Both functions and A^? have extremum 
placed at the same distance depending on x and atomic number A. This is a consequence of the 



scaling behavior (|4.15|) and equations ( |4.16|) and (|4.17 ). The extremum occures at certain r, 



' max i 



such that N%'(T max ) = -T max N%"(T max ). In Fig. |, T max is approached by varying r± at fixed x. 
Alternatively it can be reached by varying x at fixed rj_. 

The position of the maximum r max can be used as another definition of saturation scale: 
• Definition (d): 

r(r ± = 2/Qf) = r max . (4.19) 

The saturation scale estimated from definition (d) is presented in Figs. |2] and [| 

For the sake of briefness we skip plots representing the ratios /N® and 3?/A/jP. The r_|_ 
independence of the ratios is approximately reproduced. Thus the scaling with respect to x 
variable is established. Within relative error of order 20% the resulting ratios do not depend 
neither on x no on the atomic number A. This observation is consistent with the power law 
ansatz for the saturation scale ( ft.l4|) and the conclusions of the previous section. The results on 



the scaling practically do not alter when x is varied. 



5 a diff /a tot 

In this section we consider a ratio of the inclusive diffractive dissociation to the total inclusive 
production. The mass maximally produced in the diffractive process can be related to the minimal 
rapidity gap Y = lnl/a;o: 

M 2 = Q 2 (x /x - 1) . 



6 



Ne x = 10" 



Ne x = 10" 



Ne x = 10" 



Ne x = 10" 6 



0.15 
0.1 

0.05 


-0.05 
-0.1 
-0.15 



0.2 
0.15 

0.1 
0.05 


-0.05 
-0.1 



0.2 
0.15 

0.1 
0.05 


-0.05 
-0.1 



/ 

/ 








2 4 


6 8 


10 12 


Mo x 


= 10" 


-3 


/ 

1 








2 


4 


6 8 


Au x 


= 10" 


-3 


/ 

/ 








2 


4 


6 8 



rJGeV 1 ) 



0.15 
0.1 
0.05 


-0.05 
-0.1 
-0.15 



0.2 
0.15 

0.1 
0.05 


-0.05 
-0.1 



0.2 
0.15 

0.1 
0.05 


-0.05 
-0.1 



/ 






2 4 6 8 


10 12 


Mo X = 10" 


-4 


/ ^ \ 

N 

/ 

/ 








2 4 

Au x = 10" 


6 8 

-4 


/ 

/ 








2 4 


6 8 



rJGeV) 



0.15 
0.1 
0.05 


-0.05 
-0.1 



0.15 
0.1 
0.05 


-0.05 
-0.1 



0.15 
0.1 
0.05 


-0.05 
-0.1 



i \ 

I 






2 4 


6 8 10 12 


Mo X 


= 10~ 5 


r \ 
1 \ 

1 


\ 




2 


4 6 8 


Au x 


= 10~ 5 


1 \ 

1 


\ 




1 2 


3 4 5 6 



r 1 (GeV" 



0.15 
0.1 
0.05 


-0.05 
-0.1 



0.15 
0.1 
0.05 


-0.05 
-0.1 



i \ 

I 


\ 






1 


2 3 


4 5 6 7 


Mo 


X = 


io~ 6 


I \ 

I 


\ 






1 


2 


3 4 5 


Au 


X = 


10" 6 



r^GeV" 1 ) 




Figure 4: T/ie derivative functions (dashed line) and (solid line) as functions of transverse 
distance at fixed xq = 1CT 2 . 
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At x = Xq diffraction reduces to pure elastic scattering. We set Xq = 10 2 as a minimal gap 
allowed in our calculations. 

Fig. [| displays the ratio a d% ^ ja tot as a function of x for fixed values of the photon virtuality Q 2 . 
This ratio measures the value of shadowing corrections and it grows as x decreases tending to the 
unitarity bound 1/2. For heavier nuclei, diffraction dissociation is larger as a result of stronger 
nuclear shadowing. 
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Figure 5: The ratio of inclusive diffractive dissociation to total inclusive production is plotted versus x. 
The curves correspond to Au (solid line), Mo (dashed line), and Ne (dotted line) nuclei. 

It is worth to investigate the A-dependence of cr dl fl jo tot which in perturbative QCD is propor- 
tional to A 1 / 3 (times a dl ^ / a tot for a proton). We assume the ratio a dl ^ ' / a tot to have the power 
law dependence on A: 

a diff j a tot ^ A n{x) _ 

It is clear that in deep saturation regime n should vanish leading to the A-independent behavior 
(a d,l f* ja tot = 1/2). Fig. || displays the function n(x) for all nuclei in consideration. 



n 




Figure 6: The exponent n is plotted ver- 
sus x for several values of Q 2 : 
Q 2 = 80 GeV 2 (solid line), 
Q 2 = WGeV 2 (dotted line), 
and Q 2 = I GeV 2 (dashed 
line). 
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It is important to note that the results obtained are in agreement with the results of Ref. pi 
where the ratio a d% ^ ja tot was considered on a basis of the Glauber-Mueller formula. 



8 



6 Theoretical discussion. 



One of the most surprising results of our calculations is the geometrical scaling behavior which 
holds even at sufficiently short distances (large values of photon virtualities Q 2 ). At short distances 
the imaginary part of the dipole elastic amplitude has the form 

4 i 4 

N(A;r ± ,x) oc r 2 ± xG A (x, ^-)/ttR 2 a « As r]_ xG proton (x, -=-). 

The scaling means that this function is a function of only one variable r = rj_Q 2 (A; x). Such a 
nontrivial behavior should certainly result from a kind of interplay between three variables r±, x 
and A. For a proton target it was shown in Ref. p2| that the scaling holds in a wide range of 
distances even for large virtualities Q 2 ss 4/r 2 L > Q 2 (x): 

Q 2 < Q 2 < Q A JA 2 ■ (6.20) 



However, for the geometrical scaling to hold in the case of nuclear targets, an additional condition 
on the value of the saturation scale is required. Namely, the condition Q 2 S (A; x) ^> A 2 As should 
be fulfilled in order to justify the interplay between A and the rest of the variables. As was shown 
in Ref. 
follows 



23|j, in the region of low x and large photon virtualities the dipole amplitude looks as 



N(A;£ = ln(l/(r 2 ± A 2 )),y - y = ln(x /x)) « e 2 V«^(«») + ^^ , (6.21) 
where £ A = g In A. 

Saturation scale can be found from the condition N = Const(A, r±, x), which leads to the 
equation for the saturation scale £ s = \n(Q 2 / A 2 ): 



4a s Ut{y-yo) = (Ut - U) 2 ■ (6.22) 

As a result, £ sat (y = y ) = £ A while £ sat = 4 a s (y - y ) + 2£ A for £ sai > £ A . Let us expand 
Eq. ( ggjj ) in the vicinity of Q 2 « Q 2 . Introducing £ = A£ + £ sat with A£ = ln(Q 2 /Q^ and 
considering A£ <C £ sai we get 



2^a s £(y-y ) + U - £ = ~A£ + A£(^_) + O ( (f^) 2 . (6.23) 

We are interested in the region where an approximate geomentrical scaling behavior can be seen. 
So, we consider A£ « 1 and for the scaling behavior to hold we have to assume that 

A£ < £ sa t and £ A < 6at • 
The first inequality leads to Q 2 < Q\j 'A 2 while the second one gives Q 2 > A 2 A3. 



9 



In the particular model given by Eq. ( |6.21|) and at low x, the A dependence of the saturation scale 



squared is Q 2 oc Aa. For such a dependence the condition Q 2 > A 2 A3 means that x is supposed 
to be sufficiently small. Our numerical solution, however, shows rather different dependence on A 
[Q 2 oc As) and the requirement Q 2 > A 2 A3 has numerical justification only. We would like to 
recall that the above model is correct at very large photon virtualities only. A natural question 
arises how to extend this model to smaller virtualities in the range defined by Eq. ( |6.2(J| ). We will 
discuss this issue below. 

The above discussed simple model allowed us to illustrate that the theoretical expectations for 
the A-dependence of the saturation scale are quite different from the numerical calculations. A 
possible reason for such different A-dependence is in the fact that x = 10~ 2 — 10~ 5 are not small 
enough to apply this model. For rather low energies the saturation scale squared obtained from 
Eq. (|6.22|) is proportional to A^ in agreement with the numerical calculations. In fact, we are still 
far away from the low x region where the collective effects related to saturation are strong. This 
can be seen from the ratio a dl ^/a tot which is not close to the saturation limit 1/2 and moreover 
exhibits some A dependence. 

Now, let us try to analyze the diffraction dissociation in DIS using the same approach as in the 
discussed model. We consider so large initial virtualities that we can restrict ourselves to solutions 
of the linear DGLAP evolution equation. The diagram which contributes to high mass diffraction 
is shown in Fig. [7[ It is well known (see Ref. [24] for example) that this diagram can be expressed 
in the form 

N D (y-y = \nM 2 ,0 = J d?N(y - y ,£ - £') N 2 (y ,£) (6.24) 
with N(y - y , £ - f ') given for f - f ' » £ strf by 

= /^- e t+(/- 1 )(^') oc e V a «(!/-i/o) + . (6 . 25 ) 
J 2ni 



It can be seen that Eq. ( |6.24j ) is a solution to Eq. ( |1.5| ) provided the non-linear terms are 



neglected in this equation. Moreover, Eq. (|6.24j ) satisfies the correct initial conditions. The non- 



linear corrections are neglected since we wish to estimate the saturation scale which corresponds 
to sufficiently small A^ D (in semiclassical approach, for example, N D oc as). It is important to 
emphasize that the A dependence in Eq. fl6.24j ) comes from the initial conditions only. 

Assume the value of xq to be so small that N(yo,l;') is saturated at the saturation scale £' = 
t,sat{yo)- In this case, the typical value of £' in the integrand of Eq. (|6.24|) is £' ~ £ S at(yo)- Indeed, 



for £' < £ sa t(yo), N(y ,£') ~ const(y ,^') and the main contribution to the integral comes from 
the upper limit of the integration £' = Csat(vo)- in t ne region £' > C, sa t(yo) the integrant falls 
down exponentially as a function of As a result, the diffractive amplitude A^ D is approximately 
equal 



N D (y - J/,, = hiM 2 ,0 = N(y -y ,Z- Ut(yo)) N 2 (y ,Ut) ■ (6.26) 
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Since N(y ,£ sat ) = Const(y ) we have 



N D (y-y = hxM 2 ,£) oc e V a « ) (w )) -€+£..* too ) _ ( 62 7) 
Eq. ( |6.27D leads to the saturation scale given by the equation 



4a s (y - 2/ ) = CM - Ut(Vo) ■ (6.28) 

Recall that £ 8a t(yo) = 4a s y + £a (£«at(yo) ~ 60- Finally we obtain from Eq. (|6.28|) that 

C D at (y)=4a s2/ + U- (6-29) 

It can be seen that this simple approach reproduces our numerical result that the saturation scale 
for the diffractive production in DIS has the very same x dependence as the saturation scale for 
the total DIS cross section. However, it is important to note that the energy dependence of the 
saturation scale obtained from these simple theoretical estimates turns out to be quite different 
from the numerical results. Within the approximations made, £,® at (y) does not depend on y (x ) 
in accordance with the numerical calculations. 



Comparing Eqs. ( |6.22j ) and ( |6.28j ), the A-dependence of the saturation scale in total and diffractive 



productions appear to be quite different. Technically, this difference arises from the extra £' 
integration in Eq. fl6.24| ) in comparison with Eq. ( |6.21[ ). Keeping this in mind we are going to 
reconsider saturation in total cross section. To this goal, we can rewrite the dipole amplitude iV 
as an integral over the intermediate virtuality As follows from the general properties of the 
DGLAP equation this can be always done. 

Fig. H displays the DIS process in which the initial condition for the DGLAP linear evolution is 
fixed by the McLerran-Venugopalan model (see Ref. [Q]). In other words, it is assumed that at 
sufficiently low energies (y' in Fig. ||), saturation at Q 2 (A,y') ^> A 2 is reached due to strong 
gluonic fields in a nuclear target. In this case the integral over £' is dominated by £' ~ C, sa t(y') 
and the the dipole amplitude N is given by 

N(A;Z = \n(l/(rlA 2 ),y-y f = \n(x'/x)) sa e 2 V / ^ I ^M + Ut^)-* . ( 6 . 30 ) 

Since £, sa t{y') ~ £a i n the initial condition for the saturation scale, the following expression is 
obtained: 

Ut(y) = 4a s (y - y') + Ut(y') = 4a s (y - y') + U . (6.31) 



The saturation scale defined in Eq. ( |6.31| ) is proportional to A3 in a contrast to Eq. ( |6.22| ). A 



natural question arises: what is wrong with the first approach? Eq. ( |5.21 ) satisfies all correct 



initial conditions at short distances. In the integration over £' the main contribution comes from 
the upper limit of the integration £' w t; sat . The integral over in the saturation region looks as 
follows: 
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y* 



yo^sat(y(P 



N(y-y .?- Uv )) 



n (y . ? sa ,(y)) 



Figure 7: T/ie diffractive production of large mass (M) (y — y$ = \nM 2 ) in DIS. 




N(y-y 



n ( y , $Jy)) 



Figure 8: The total cross section for DIS with initial conditions from the McLerran- Venugopalan model. 
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For £' > £ sat the integral is d£' e"^'. This simple form of the integral is valid in the region 
Q 2 > Qg/A 2 . This contribution is small but the estimates in the beginning of this section show 
that at sufficiently short distances defined by Eq. (|6.22 ) even such a small contribution may reach 



so large values that we have a dense system of partons. 
Actually, Eq. is 



N(A; f = hi(l/(r| A 2 ),y-y = ln(x /x)) =r\A xG(x, \)/vR\ (6.32) 

r < 



in usual notations with xG being in double log approximation of pQCD. Eq. ( |6.32| ) is correct at 
least at very short distances. 

Within the same notations Eq. ( p. 30 ) looks differently 



N(A; £ = ln(l/(rl A 2 )),y - y = \n(x /x)) =r\A xG(x, * )/t:R\ . (6.33) 

r±Qt{A;y') 



As we discussed (see Ref. for more details) the geometrical scaling behavior of Eq. ( |6.33D 
is preserved till r\ ps 4A 2 /(As Q*(y')). For shorter distances we expect the regime discussed in 
the beginning of this section to take place with Eq. fl6.32| ) being correct. This should happen at 
Q 2 > Q sat / '(As A 2 ) and at x — > 0. Unfortunately all experimentally accecible values of x are not 
small enough for Eq. ( p.21] ) to be seen in a nearest future. 

To evaluate a typical scale of distances involved we recall that the first RHIC data are likely to 
reveal Q 2 {Go\d, x = 10~ 2 ) ~ 2 GeV 2 |p5| . Therefore, the regime with Q 2 oc A3 is expected to 



start at r\ < 0.01 GeV 2 = 4 10 4 fm 2 . As can be seen in Fig. [[| at so short distances and at 
x = 10~ 2 — 10~ 6 the diffraction amplitude is far away from the saturation region. 



7 Conclusions 



The non-linear evolution equation ( |1.5|) is solved numerically for six real nuclear targets Ne2o, 
Ca4o, Zn 70 , Mo 100 , NdisQ, and Au 197 . The obtained solutions display the very same pattern as 
in the case of proton target 



16fl . These solutions are used to study the A-dependence of single 



diffractive dissociation. 



The saturation scale Qj A is estimated. The fit to the parameterization Qj(A, x, Xq) ~ A p x~ 
reveals powers p and A. The results on A coincide with the corresponding power obtained for 
the proton case [16] while the A-dependence is the same as found for the saturation scale Q s ,a 
- saturation scale obtained in total production JT2" |. The function Qf A is found to be almost 
independent on x . 



The geometrical scaling with respect to x is well established for all nuclei considered. The scaling 
holds within a few percent accuracy and in the whole kinematic region investigated. As a con- 
sequence, inclusive diffractive production off nuclei is predicted to possess the scaling both with 
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respect to x and A. Like for a proton target the scaling phenomena with respect to Xq set in at 
x-Cio but is violated at x ~ x . 

The A and Xq dependences of the determined saturation scale can be given a theoretical expla- 
nation. However, the numerically found x dependence of the scale is weaker compared to the 
theoretical estimates. 

The ratio a diff /a tot was examined on a basis of the solutions obtained. At x ~ 10 4 a significant 
shadowing is expected of the order 25% and it is larger for heavy nuclei. The fact that this ratio 
turns out to be very close to the estimates based on the Glauber-Mueller formula shows that 



the latter can be used as a simple approach for first estimates of possible collective effects in the 
saturation region of high parton density QCD. 
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